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The complete elastic tensor of YM2B2C was determined by application of the resonant ultrasound 
spectroscopy technique to a single-crystal sample. Elastic constants were found to be in good 
agreement with partial results obtained from 'time-of-flight' measurements performed on samples 
cut from the same ingot. From the measured constants, the bulk modulus and Debye temperature 
are calculated. 



The discovery of superconductivity in .the borocar- 
bide quaternary intermetallic compoundsu has gener- 
ated a great deal of interest due to the observed in- 
terplay between superconductivity and magnetic order. 
Compounds in this family have formulae of the form 
i?Ni2B2C, where choice of the rare earth ion, R, strongly 
affects the physical properties, which range from su- 
perconducting with no magnetic ordering (R—Y) to 
non-superconducting with a Neel temperature of 19 K 
(i?=Gd)E 

Recent heat transport measurements on non-magnetic 
borocarbides suggest that the superconductivity is veriz. 
anomalous, with a .gap that is either highly anisotropico 
or has point nodes. □ To elucidate the nature of the order 
parameter, a powerful tool is ultrasound attenuation (see 
for example Ref. ||). As a first step towards measuring 
the attenuation of sound it is important to establish the 
full elastic tensor of the material. 

The 2 nd order elastic constants relate stresses placed 
on a material along various directions to the resulting 
deformations of the crystal lattice. As such, knowledge 
of a compound's elastic constants can be useful when in- 
vestigating a number of solid state properties, such as in- 
teratomic potentials, equations of state, phonon spectra, 
the Debye temperature, and the bulk modulus. In the 
most general case, a highly anisotropic solid can have 21 
independent second-order elastic constants, but the crys- 
tal structure for all compounds in the iJN^P^C group is 
body-centered-tetragonal (space group 14/ 771mm), con- 
sisting of R-C planes separated by N12B2 layers stacked 
along the c axis,El reducing the number of independent 
constants to 6.13 In Voigt notation, these constants are 
denoted Cn, C33, C13, C%2, C44, and C§§; all other 
tensor entries can be retrieved from these six. 

In this article, the elastic constants of YM2B2C are 
determined at room temperature via two methods: res- 
onant ultrasound spectroscopy and ultrasonic lime-of- 
flight. Resonant ultrasound spectroscopy (RUS)Q allows 
the determination of the complete elastic tensor of a ma- 
terial in one measurement, and was used as the primary 
experimental technique of this study. As a complement to 
the RUS measurements, elastic constants were also deter- 
mined by measuring sound velocities using the more con- 
ventional 'time-of-flight' method. Note that this study is 
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FIG. 1: Schematic of a typical RUS experimental set-up (from 
Ref. |2l). 



not simply a confirmation of existing measurements: it 
adds a missing elastic constant and provides a very ac- 
curate second measure to secure time-of-flight results in 
slight disagreement. 

The RUS method utilizes the normal modes of vibra- 
tion of a solid sample of known density, geometry and 
crystal structure to determine the elastic tensor. To per- 
form an RUS measurement, the normal mode frequency 
spectrum of the sample in question is first obtained. This 
is done by holding the sample between two piezo-electric 
transducers, generally by its corners in order to reduce 
damping of the vibrational modes. One transducer is 
driven at a frequency which is swept across the region 
of the spectrum containing the normal modes, thus driv- 
ing the sample resonances, while the other picks up the 
resulting amplitudes of crystal vibration (see Fig. ^) . 

Once a list of normal mode frequencies has been com- 
piled, the elastic constants can be calculated. While it is 
possible to directly compute the normal mode frequen- 
cies of a sample, given the dimensions, mass and elas- 
tic constants, there is no known analytic method to do 
the reversa, so a computerized fitting algorithm has been 
developed!] to numerically complete the calculation. This 
algorithm changes the elastic constants, calculates result- 
ing normal mode frequencies, compares these calculated 
frequencies to the experimentally measured list and then 
iteratively repeats the process until the best fit to the 
data is obtained. 

In the time-of-flight method, short ultrasound pulses 
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are propagated through a sample carefully prepared with 
two flat, parallel faces, and later detected in order to ob- 
tain the sound velocity in the direction of pulse prop- 
agation. The diagonal entries of the elastic tensor can 
be easily calculated using the velocities measured along 
various principal crystal directions. The off-diagonal con- 
stants, however, require measurements of sound velocities 
in non-principal directions, which can be difficult. The 
equation (Cu + Ci 2 + 2C m )/2 = p v 2 [ 110 / L ], where p 
is the mass density of the material, and «[ho/l] is the 
velocity of sound propagating in the [110] direction with 
a longitudinal polarization, is an example of how a sound 
velocitv-js related to the elastic constants for a tetragonal 
crystalB In this notation scheme, a transverse polariza- 
tion, for example in the [010] direction, is indicated by 
T010. 

With the complete elastic tensor, the bulk modulus of 
a tetragonal crystal, Bo = — V(^), is easily obtained 
via the equation: 



C^(C\\ + C12) — 2C 



C\\ + C12 + 2C; 
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4C 
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(1) 



Also, an estimate of the Debye temperature can be 
made. For axrystal with tetragonal structure the follow- 
ing equations is employed: 



id 



I: (9NY 



-1/3 



p- l / 2 J- 1 ' 2 , 



(2) 



wherein h is Planck's constant, Ub is Boltzmann's con- 
stant, N is the number of atoms per unit cell, V is the 
volume of the unit cell, p is the mass density of the ma- 
terial, and J is the harmonic series expansion of an inte- 
gration over sound velocities in all directions, written in 
terms of the elastic constants (Eqs. 32 and 33 in Ref. ||). 

A large high-quality single crystal specimen of 
YNi2 n B2C was grown using the floating zone method. 
From this specimen a sample was cut using spark erosion 
and then polished into the shape of a rectangular paral- 
lelepiped, with sides aligned along the principal crystal 
axes as determined by Laue x-ray backscattering mea- 
surements. 

Dimensions of the polished sample were measured to 
be 0.0633(8) cm x 0.0873(9) cm x 0.238(2) cm. The 
sample mass was 0.00832(3) g, yielding a density of 
6.3(1) g/cm 3 . Note that this is somewhat different 
from the density of 6.09(1) g/cm 3 for YNi2 n B2C calcu- 
lated using latti^canatants. obtained from x-ray/neutron 
measurementsMBByy 

A number of other measurements do indicate good 
sample quality. Specifically, a residual resistivity ratio 
(RRR) of 14.2, a superconducting width of 0.3 K, and 
a T c of 14.7 K, were observed via resistivity measure- 
ments on our specimen,, wh^h Cf' T 7?Pfj' T T favourably with 

published vaiuesEmBemSm 

The experimental apparatus used to detect normal 
mode frequencies was a DRS Modulus I Resonant Ultra- 
sound Spectrometer, with lead zirconate titanate trans- 



ducers. This system has previously been used to success- 
fully determine the elastic constants of Sr2Ru04.Eil 

Once compiled, the list of normal mode frequencies was 
fitted, and appropriate elastic constants obtained. The 
measured dimensions of the sample were used as initial 
dimensions for the program input, but as the measuring 
errors in these values were the largest sources of error in 
this study, the dimensions were left as free parameters 
in the fit. Thus, 9 free parameters - 6 elastic constants 
and 3 dimensions - were fitted to the measured frequency 
data, while quantities such as the mass and density of the 
sample were held constant. 

Two other samples were cut from the same ingot as 
the RUS sample and polished to facilitate time-of-flight 
measurements in the [100] and [110] directions. Sam- 
ple thicknesses were measured to be 0.319(2) cm and 
0.227(2) cm for the samples oriented along the [100] and 
[110] directions, respectively. In all cases, transducers 
were bonded to samples using phenyl salicylate and sil- 
icone for room and low temperature measurements re- 
spectively. Time-of-flight measurements were made at 
room temperature (300 K) and low temperature (2 K) 
using 20 and 30 MHz LiNbC>3 transducers and a home- 
built spectrometer £3 which utilizes the phase comparison 
technique at frequencies between 20 and 500 Mhz. This 
apparatus has previously been used irua successful study 
of ultrasonic attenuation in Si'2Ru04.El 

Using the RUS technique, a total of 68 normal modes 
were measured between 0.5 MHz and 4.5 MHz. The fre- 
quencies of these measured modes, along with the fre- 
quencies of the calculated modes from the best fitting 
run and the percent difference between them, are shown 
in Table [| The corresponding frequency spectrum is 
shown in Fig. [2[ with a representative portion of the 
spectrum inset in more detail. In the full-range view, 
a large amount of noise was measured between approxi- 
mately 0.7 MHz and 1.1 MHz, and is attributed to fre- 
quency characteristics of the measuring apparatus itself. 
A prototype RUS apparatus, with PVDF (polyvinylidene 
fluoride) film transducers JH3 was used to verify that no 
normal mode frequencies of the sample lay in this range. 
In addition, the fits indicate no missing modes in this 
region. 



The best fit to the data had a root-mean-square error 
between measured and calculated frequencies of 0.38%, 
indicative of a good fit. The fitted dimensions were 
0.0636(6) cm x 0.0878(8) cm x 0.236(2) cm, which 
agree with measurements, also providing confirmation of 
a good fit. Note that only two modes were not detected, 
but were determined to be existent through the fitting 
procedure. Also, the second mode was not weighted in 
the fit, a common practice among the first 1-2 modes in 
RUS measurements. U 

The elastic constants determined via the best fit are 
listed in Table |h[ Error values are obtained from the 
largest of either an estimation of the quality of the fit, as 
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TABLE I: Comparison between experimentally measured normal mode frequencies and normal mode frequencies generated by 
fit. 
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4.514518 


4.520470 


0.13 



"not weighted in the fit, as discussed in the text. 

'missing modes, not detected by the apparatus as explained in the text. 



generated by the fitting program, or from the variation 
in Cij values over the range of densities allowed by the 
error in the measured dimensions. Specifically, the errors 
in C44 and all three fitted dimensions come from the 
error in density, while the rest come from quality-of-fit 
estimates. 

Also included in Table || arc the clastic constants de- 
termined from three different sources using the time-of- 
flight technique: room temperature measurements from 
this study, low-temperature measurements done by our 



groups and low-temperature measurements by Isida et 
al.c3 The elastic constants from our room temperature 
and low-temperature time-of-flight measurements j-svere 
calculated from velocity data as discussed above.El An 
example of the raw time-of-flight data used to obtain ve- 
locities is shown in Fig. ||. The time delay between each 
echo is equivalent to the total travel time of an ultra- 
sound pulse going through the sample and back; when 
taken with the thickness of the sample this yields the 
velocity of sound for that particular direction of propa- 
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1 2 3 4 5 

Frequency (MHz) 

FIG. 2: The frequency spectrum measured for YN12B2C. Inset: a representative portion of this spectrum is expanded, covering 
a region from 2.85 MHz to 2.96 MHz in detail. 




20 30 
time (us) 

FIG. 3: Example of raw time-of-flight data, in the case of a 
low-temperature (2 K) measurement of V[ 110 /tooi] ■ 



gation/polarization combination. 

The measured and calculated sound velocities with the 
appropriate polarization and direction are listed in Ta- 
ble Note that whenever the density of YNI2B2C was 
required in a calculation, p — 6.3(1) g/cm 3 was used for 
measurements performed on our specimen, and p = 6.05 
g/cm 3 was used for the measurements by Isida et al. 

Upon inspection of the data shown in Tables |0] and 111 , 
it is immediately apparent that while the values obtained 
via both the RUS and time-of-flight methods of our grooip 
agree within error, the values published by Isida et alx3 : 
with the exception of C^e, are between 16% and 36% 
lower than the RUS results and both sets of time-of-flight 
results from our group. Reasons for this disagreement are 



not clear, but agreement of RUS and time-of-flight mea- 
surements on samples cut from the same specimen may 
suggest the significance of differences in sample compo- 
sition between the two groups. 

The RUS values of the elastic constants give a bulk 
modulus for YNi 2 B 2 C of 1.8(1) xlO 12 dynes/cm 2 . This 
is in the range of the isothermal bulk modulus values of 
1.2(1) xlO 12 dynes/cm 2 and 2.00xl0 12 dynes/cm 2 pp 
ously measured for this material by x-ray diffractionJH. 
which gives us added confidence in our results. Further- 
more, the Debye temperature was calculated via Eq. |^, 
yielding a value of 525(10) K, which falls within the wide 
range of published values of the Debye temperature for 
YNI2B2C: 340(20) KO-Jrom a resistivity measurement, 
and 415 kJ13 489(5) KB and 537 Kll3 from specific heat 
measurements . 

In conclusion, all six independent elastic constants of 
YNI2B2C have been determined using the resonant ultra- 
sound spectroscopy method. Good agreement was found 
between these results and the elastic constants measured 
using the time-of-flight technique on a sample from the 
same ingot. A calculation of the bulk modulus and an es- 
timation of the Debye temperature of YM2B2C were per- 
formed with the elastic constants obtained in this study, 
and the results were found to be consistent with the range 
of values found in the literature for these quantities. 



Acknowledgments 

This work was supported by the Canadian Institute 
for Advanced Research and funded by NSERC. J. P. also 
acknowledges the support of the Walter C. Sumner Foun- 
dation. 



5 





TABLE II: Elastic constants of single crystal YM2B2C, 


in units of 10 12 


dynes/cm 2 . 
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14.2 


"from Ref. list 












'calculated from the quoted values of Cn and (On - Ci2)/2 in Ref. |l8| 












TABLE III: Sound velocities in YNi 2 B 2 C, in units of km/s. 






Method 


V[100/L] 


V [001/L] V[no/T001] , V[ioo/T001] V[ioo/T010] 


V [110/T1101 




Temperature (K) 


RUS 


6.8(2) 


6.4(2) 3.20(5) 4.74(8) 


3.3(3) 


7.6(2) 


300 


time-of-flight 


6.72(6) 


3.26(2) 4.76(4) 


3.36(2) 


7.66(4) 


300 


time-of-flight 


6.8(3) 


3.3(2) 4.5(6) 


3.3(2) 


7.4(4) 


2 


time-of-flight" 


6.0(3) 


5.9(3) 3.0(1) 4.7(2) 


3.2(2) 


6.9(7) 


14.2 



from Ref. |l8[ 
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